REPORT  DOCUMENTATION  PAGE 


Form  Approved  0MB  NO.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments 
regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggesstions  for  reducing  this  burden,  to  Washington 
Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington  VA,  22202-4302. 
Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  oenalty  for  failing  to  comply  with  a  collection 
of  information  if  it  does  not  display  a  currently  valid  0MB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 

2.  REPORT  TYPE 

3.  DATES  COVERED  (From  -  To) 

06-07-2014 

Final  Report 

1-Apr-2012-31-Mar-2013 

4.  TITLE  AND  SUBTITLE 

Role  of  Nanometer-Scale  Mechanical  Responses  of  Hybrid 
Ceramic-Based  Materials  on  Toughening 


5a.  CONTRACT  NUMBER 

W911NF-12-1-0153 


5b.  GRANT  NUMBER 


6.  AUTHORS 
Susan  B.  Sinnott 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMES  AND  ADDRESSES 

University  of  Florida  -  Gainesville 
219  Grinter  Hall 
PO  Box  115500 

Gainesville,  FL  32611  -5500 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS 
(ES) 

U.S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-221 1 


12.  DISTRIBUTION  AVAILIBILITY  STATEMENT 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


10.  SPONSOR/MONITORS  ACRONYM(S) 
ARO 


1 1 .  SPONSOR/MONITORS  REPORT 
NUMBER(S) 

62034-MS-DRP.l 


Approved  for  Public  Release;  Distribution  Unlimited 


13.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  contrued  as  an  official  Department 
of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


14.  ABSTRACT 

We  used  classical  molecular  dynamics  simulations  to  identify  the  nanometer-scale  mechanisms  by  which  a 
ceramic-based  hybrid  composite  of  PMMA  and  A1203  responds  to  mechanical  deformation.  The  influence  of 
factors  such  as  the  arrangement  of  the  phases  in  lamellar  and  brick-and-mortar  structures,  atomic-scale  roughness 
of  the  A1203  phases,  and  attachment  of  molecular  grafts  to  the  A1203  were  explored.  We  also  used  first-principles, 
density  functional  theory  calculations  to  examine  the  way  in  doping  of  yttrium  aluminum  garnet  (YAG)  influenced 


15.  SUBJECT  TERMS 

Hybrid  ceramic-polymer  materials,  YAG,  molecular  dynamics  simulations,  density  functional  theory  calculations,  mechanical 
properties,  optical  properties 


Susan  Sinnott 


19b.  TELEPHONE  NUMBER 
352-846-3778 


Standard  Form  298  (Rev  8/98) 
Prescribed  by  ANSI  Std.  Z39. 1 8 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

15.  NUMBER 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF  PAGES 

uu 

UU 

uu 

UU 

Report  Title 

Role  of  Nanometer-Scale  Mechanical  Responses  of  Hybrid  Ceramic-Based  Materials  on  Toughening 

ABSTRACT 

We  used  classical  molecular  dynamics  simulations  to  identify  the  nanometer-scale  mechanisms  by  which  a  ceramic-based  hybrid  composite 
of  PMMA  and  A1203  responds  to  mechanical  deformation.  The  influence  of  factors  such  as  the  arrangement  of  the  phases  in  lamellar  and 
brick-and-mortar  structures,  atomic-scale  roughness  of  the  A1203  phases,  and  attachment  of  molecular  grafts  to  the  A1203  were  explored. 
We  also  used  first-principles,  density  functional  theory  calculations  to  examine  the  way  in  doping  of  yttrium  aluminum  garnet  (YAG) 
influenced  its  optical  properties.  In  particular,  the  identity  of  the  dopants,  their  location  within  YAG  unit  cell,  and  their  concentrations  were 
considered  with  guidance  from  experimental  data  from  Lawrence  Livermore  National  Laboratory. 


Enter  List  of  papers  submitted  or  published  that  acknowledge  ARO  support  from  the  start  of 
the  project  to  the  date  of  this  printing.  List  the  papers,  including  journal  references,  in  the 
following  categories: 

(a)  Papers  published  in  peer-reviewed  journals  (N/A  for  none) 


Received  Paper 


TOTAL: 


Number  of  Papers  published  in  peer-reviewed  journals: 


(b)  Papers  published  in  non-peer-reviewed  journals  (N/A  for  none) 


Received  Paper 


TOTAL: 


Number  of  Papers  published  in  non  peer-reviewed  journals: 


(c)  Presentations 

“Investigation  of  A1203  Nanostructure  Surfaces  Using  Charge  Optimized  Many  Body  Potentials”,  D.  Yilmaz,  T.  Liang,  S.  Phillpot,  S. 
Sinnott,  AVS  59th  International  Symposium  and  Exhibition  in  Tampa,  Florida,  October  29-November  2,  2012. 


Number  of  Presentations:  1 ,00 

Non  Peer-Reviewed  Conference  Proceeding  publications  (other  than  abstracts): 


Received  Paper 


TOTAL: 


Number  of  Non  Peer-Reviewed  Conference  Proceeding  publications  (other  than  abstracts): 

Peer-Reviewed  Conference  Proceeding  publications  (other  than  abstracts): 


Received  Paper 


TOTAL: 


Number  of  Peer-Reviewed  Conference  Proceeding  publications  (other  than  abstracts): 

(d)  Manuscripts 


Received  Paper 


TOTAL: 


Number  of  Manuscripts: 

Books 


Received  Book 


TOTAL: 


Received 


Book  Chapter 


TOTAL: 


Patents  Submitted 


Patents  Awarded 


Awards 

Fellow  of  the  Materials  Research  Society,  2012 
Fellow  of  the  American  Physical  Society,  2013 


Graduate  Students 

NAME 

PERCENT  SUPPORTED  Discipline 

Kai  Yao 

0.33 

Kamal  Choudhary 

0.66 

FTE  Equivalent: 

0.99 

Total  Number: 

2 

Names  of  Post  Doctorates 


NAME 

PERCENT  SUPPORTED 

Dundar  Yilmaz 

1.00 

FTE  Equivalent: 

1.00 

Total  Number: 

1 

Names  of  Faculty  Supported 


NAME 

PERCENT  SUPPORTED  National  Academy  Member 

Susan  Sinnott 

0.14 

FTE  Equivalent: 

0.14 

Total  Number: 

1 

Names  of  Under  Graduate  students  supported 


NAME 

PERCENT  SUPPORTED 

FTE  Equivalent: 

Total  Number: 

Student  Metrics 

This  section  only  applies  to  graduating  undergraduates  supported  by  this  agreement  in  this  reporting  period 

The  number  of  undergraduates  funded  by  this  agreement  who  graduated  during  this  period: . 0.00 

The  number  of  undergraduates  funded  by  this  agreement  who  graduated  during  this  period  with  a  degree  in 

science,  mathematics,  engineering,  or  technology  fields: . 0.00 

The  number  of  undergraduates  funded  by  your  agreement  who  graduated  during  this  period  and  will  continue 

to  pursue  a  graduate  or  Ph.D.  degree  in  science,  mathematics,  engineering,  or  technology  fields: . 0.00 

Number  of  graduating  undergraduates  who  achieved  a  3.5  GPA  to  4.0  (4.0  max  scale): . 0.00 

Number  of  graduating  undergraduates  funded  by  a  DoD  funded  Center  of  Excellence  grant  for 

Education,  Research  and  Engineering: . o.OO 

The  number  of  undergraduates  funded  by  your  agreement  who  graduated  during  this  period  and  intend  to  work 

for  the  Department  of  Defense . 0.00 

The  number  of  undergraduates  funded  by  your  agreement  who  graduated  during  this  period  and  will  receive 

scholarships  or  fellowships  for  further  studies  in  science,  mathematics,  engineering  or  technology  fields: . 0.00 

Names  of  Personnel  receiving  masters  degrees 

NAME 
Kai  Yao 

Total  Number:  1 

Names  of  personnel  receiving  PHDs 

NAME 

Total  Number: 

Names  of  other  research  staff 

NAME  PERCENT  SUPPORTED 

FTE  Equivalent: 

Total  Number: 


Sub  Contractors  (DD882) 


Inventions  (DD882) 


Scientific  Progress 


Summary  of  the  most  important  results 

In  the  case  of  the  hybrid  material  project,  the  most  important  results  include: 

•  Development  of  computational  method  to  enable  the  building  of  hybrid  polymer-ceramic  structures  without  overlapping  atoms. 

•  Optimization  of  charge  optimized  many-body  potential  for  AI203  to  enable  reactive,  atomic-scale  modeling  in  simulations. 

•  Examined  laminar  and  brick-and-mortar  hybrid  material  structures  subjected  to  compression  and  tension;  documented  atomic- 
scale  responses  and  showed  how  the  responses  varied  with  the  volume  fraction  of  polymer. 

•  Introduced  atomic-scale  roughness  and  molecular  grafts  to  the  AI203  to  see  how  these  factors  changed  the  mechanical 
responses.  We  have  not  yet  been  able  to  draw  firm  conclusions  regarding  their  effects. 

In  the  case  of  the  doped  YAG  project,  the  most  important  results  include: 

•  When  Gd  dopants  are  added,  its  f-electrons  modify  the  optical  properties  of  the  YAG  in  the  IR  region.  In  contrast.  Sc  and  Ga 
are  not  predicted  to  substantially  modify  the  optical  properties  of  YAG  in  the  IR  region. 

•  Doping  with  Gd,  Sc,  and  Ga  individually  or  in  combination  with  one  another  is  predicted  to  increases  the  opaque  nature  of  the 
YAG  in  all  cases,  but  less  in  the  IR  region  than  at  higher  wavelengths. 

Technology  Transfer 


Final  Progress  Report 

Role  of  Nanometer-Scale  Mechanical  Responses  of  Hybrid  Ceramic-Based  Materials  on 

Toughening 

Susan  Sinnott 

Department  of  Materials  Science  and  Engineering,  University  of  Florida,  Gainesville,  FL  32611- 

6400 


Statement  of  the  problem  studied 

We  used  classical  molecular  dynamics  simulations  to  identify  the  nanometer-scale 
mechanisms  by  which  a  ceramic -based  hybrid  composite  of  PMMA  and  AI2O3  responds  to 
mechanical  deformation.  The  influence  of  factors  such  as  the  arrangement  of  the  phases  in 
lamellar  and  brick-and-mortar  structures,  atomic-scale  roughness  of  the  AI2O3  phases,  and 
attachment  of  molecular  grafts  to  the  AI2O3  were  explored. 

We  also  used  first-principles,  density  functional  theory  calculations  to  examine  the  way  in 
doping  of  yttrium  aluminum  garnet  (YAG)  influenced  its  optical  properties.  In  particular,  the 
identity  of  the  dopants,  their  location  within  YAG  unit  cell,  and  their  concentrations  were 
considered  with  guidance  from  experimental  data  from  Lawrence  Livermore  National 
Laboratory. 

Summary  of  the  most  important  results 

In  the  case  of  the  hybrid  material  project,  the  most  important  results  include: 

•  Development  of  computational  method  to  enable  the  building  of  hybrid  polymer-ceramic 
structures  without  overlapping  atoms. 

•  Optimization  of  charge  optimized  many-body  potential  for  A1203  to  enable  reactive, 
atomic-scale  modeling  in  simulations. 

•  Examined  laminar  and  brick-and-mortar  hybrid  material  structures  subjected  to 
compression  and  tension;  documented  atomic-scale  responses  and  showed  how  the 
responses  varied  with  the  volume  fraction  of  polymer. 

•  Introduced  atomic-scale  roughness  and  molecular  grafts  to  the  AI2O3  to  see  how  these 
factors  changed  the  mechanical  responses.  We  have  not  yet  been  able  to  draw  firm 
conclusions  regarding  their  effects. 

In  the  case  of  the  doped  YAG  project,  the  most  important  results  include: 

•  When  Gd  dopants  are  added,  its  f-electrons  modify  the  optical  properties  of  the  YAG  in 
the  IR  region.  In  contrast.  Sc  and  Ga  are  not  predicted  to  substantially  modify  the  optical 
properties  of  YAG  in  the  IR  region. 

•  Doping  with  Gd,  Sc,  and  Ga  individually  or  in  combination  with  one  another  is  predicted 
to  increases  the  opaque  nature  of  the  YAG  in  all  cases,  but  less  in  the  IR  region  than  at 
higher  wavelengths. 


Role  of  Nanometer-Scale  Mechanical  Responses  of 
Hybrid  Ceramic-Based  Materials  on  Toughening 


Prof.  Susan  B.  Sinnott 
University  of  Florida 


Final  Report 
July  6,  2014 


UF 

UNIVERSITY  of 

FLORIDA 


For  Official  Use  Only  -  Not  for  Public  Release 


UF 

UNIVERSITY  of 

FLORIDA 


Predicting  Optical  Properties  of  Doped  YAG 


For  Official  Use  Only  -  Not  for  Public  Release 


2 


Yttrium  Aluminum  Garnet 
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Y3AI2AI3O12  unit  cell  comprised  of  160 
atoms 

•  24  Y— dodecahedral 
coordination 

•  16  Al— octahedral 
coordination 

•  24  Al— tetrahedral 
coordination 

•  96  0 

Structure  with  space  group  laScf 
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Garnet  Systems 
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Computational  Methods 


Density  functional  theory  using  Vienna  Ab  initio  Simulation  Package  (VASP) 

•  PBE  pseudopotentials 

•  Single  k-point  (gamma)  calculations 

•  Cut-off  energy  =  400  eV 

•  Energy  convergence  =  10’®  eV 

•  Force  convergence  =  lO’^  eV/atom 
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Fig.  (a):  Si(0)  ~  3.40  Above:  £i(0)  ~  3.72 

Fig.  (b):  £i(0)  ~  3.54 


On  left:  Non-doped  YAG,  Xu,  Y.N.,  Chen,  Y,  et  al.,  Phys.  Rev.  B  65,  235105  (2002). 
On  right:  Results  from  this  work. 
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Electronic  structure  of  pure  YAG 
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Partial  density  of  states  (PDOS),  band-structure  of  pure  YAG 
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Introduction  of  Gd  in  YAG:  Active  in  IR  region 


Gd3Sc2Al30i2: 

Gd-f  electrons  are  contributing  extra  energy  states  in  the  infrared  region  of  YAG, 
may  be  responsible  for  infrared  activity. 
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Relation  between  heat  of  vaporization  and  energy  gap  of 
lanthanide  based  organic  and  inorganic  compounds 


•  Energy  gap  of  lanthanide  doped  YAG, 
lanthanide-TCNQ  (=  7,  7,  8,  8-tetracyano-p- 
quinodimethane)  [1]  complexes  and 
lanthanide  oxides. 

•  HOMO-LUMO  gaps  for  all  the  lanthanides 
TCNQ  complexes  were  not  available. 

•  Energy  and  HOMO-LUMO  gaps  follow 
specific  patterns. 


[1]  Zhang  et  al.,  Systematic  study  on  electrochemical 
properties  of  a  series  of  TCNQ  lanthanide  complexes. 
Journal  of  Organometallic  Chemistry  695  (2010)  1493-1498. 
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Relation  between  heat  of  vaporization  and  energy  gap  of 
lanthanide  based  organic  and  inorganic  compounds 
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Objectives 
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>  Use  classical  molecular  dynamics  simulations  to  elucidate  the 
mechanical  responses  of  PMMA-AI2O3  hybrid  structures 

>  Determine  how  the  laminar  and  brick-and-mortar  structures 
compare  to  one  another 

>  Understand  the  effects  of  degree  of  polymerization  of  the  PMMAon 
the  mechanical  properties  of  the  hybrid  material. 
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Computational  approach 
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Classical  molecular  dynamics  simulations 
Atomistic  treatment 

Forces  on  the  atoms  calculated  with  the  third-generation  charge  optimized 
many-body  potential 

Focus  on  mechanical  properties  of  AI2O3-PMMA  hybrid  materials  as  a  function 
of  hybrid  structure,  PMMA  chain  length,  and  system  size 

Currently  running:  influence  of  surface  roughness 
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Stacking  fault  energy  predicted  for  Al 
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<1 01  >  Angstrom 


<1 01  >  Angstrom 


Stacking  fault  map  predicted  for  Al  using  (left)  the  embedded  atom 
method  (EAM)  potential  and  (right)  the  charge  optimized  many-body 
(COMB)  potential. 

The  symmetry  in  the  contours  represents  the  periodicity  in  the 
system. 
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Stacking  fault  and  dislocation  splitting 
behavior  of  Al 
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Stacking  fault  3D-map  using  COMB  that  successfully  predicts  the  (121) 
direction  would  be  preferred  over(l0l)for  dislocation  propagation. 


Dislocation  splitting  of  an  edge  dislocation  within  a  (111)  plane  of  bulk  Al  is 
illustrated. 
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Mechanical  deformation  of  polycrystalline  Al  predicted 
in  MD  simulations  with  COMB  potentials 


Left:  A  13.8%  strained  Al  polycrystal 
following  constant  strain.  Red  atoms 
represent  disordered  atoms  (non-12 
coordinates),  light  blue  atoms  are  in 
an  hop  environment.  Atoms  in  an  foe 
environment  are  not  shown  for  clarity. 


(a)  (b) 


(a)  Extended  dislocations  within  the  (111) 
plane  are  indicated  by  the  black  arrows. 

(b)  The  yellow  square  highlights  twinned 
region. 
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COMB  predictions  for  AI2O3 
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Predicted  charge-spilling  across  the  Al203(0001)/AI(100)  interface 
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Uniaxial  compression  and  tension  tests  of  pure  PMMA 


PMMA  COMPRESSION 


PMMA  TENSILE 


3  4 

Strain  (%) 


Effect  of  chain  length  investigated 

Response  to  the  compression  depends  on  chain  length 

Response  to  tension  highly  non-linear 
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In-situ  brick-and-mortar  building 
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Brick  and  mortar  building  starts 
with  defining  regions  of  alumina 
(bricks)  and  PMMA  (mortar)  in 
the  simulation  box. 

First  PMMA  chains  are  built  using 
three  constraints: 

•  Angle  between  two  segments  of 
chain  should  be  greater  than 
120°. 

•  Distance  between  any  beads 
should  be  greater  than 
approximate  size  of  the  monomer 
(~0.5  nm). 

•  Distance  between  any  PMMA 
atom  should  be  larger  than  0.08 
nm  to  prevent  over-coordination 

Random  tacticity  used.  Torsion 
angle  between  monomer  and 
chain  is  set  to  avoid  overlaps 
between  polymer  atoms  in 
different  chains. 


Brick  Region 
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PMMA  Region 
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System  configuration 

deformation  axis 


Brick-and-mortar  structures  built  in  three  different  sizes  along  the  direction  of 
the  deformation  axis  with  brick  dimensions  7.8  nm,  11.7  nm,  15.6  nm  for  small, 
medium  and  large  structures,  respectively.  Other  dimensions  are  same  at  6.8  nm 
and  3.9  nm,  respectively.  There  is  a  1.2  nm  gap  between  bricks  in  all  three 
dimensions.  Periodic  boundary  conditions  are  applied  on  three  axis.  Average 
PMMA  chain  length  is  50  monomers  per  chain. 
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Illustration  of  loading  accommodation  in  the  brick-and- 

mortar  compression  test 
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Illustration  of  the  brick-and-mortar  compression  test 


ur 

UNIVERSITY  of 

FLORIDA 


Laminar  versus  brick-and-mortar  structures 


System  Dimensions: 

Bri  ck  -and  -m  o  rtar 

Small  18.4  x  8.0  x  10.4  (nm) 

Medium  26.2  x  8.0  x  10.4  (nm) 

Large  34.0  x  8.0  x  10.4  (nm) 


Lam  inar 

7.8  X  8.0  X  12.0  (nm) 
11.7  X  8.0  X  12.0  (nm) 
15,6  X  8.0  X  12.0  (nm) 


1.2  nm  gap  between  bricks,  periodic  boundary  conditions 
The  average  PMMA  chain  length  is  50  monomers  per  chain 
The  simulations  are  performed  using  the  third-generation 
Charge  Optimized  Many-Body  (COMBS)  potentials 


Initial  ^  "/<> 

Compression  of  laminar  system 
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Pressure  (MPa) 
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Uniaxial  compression  test  predictions 


Brick  &  Mortar  Compression  Test 


Laminar  Compression  Test 
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Region  1 :  Linear  region 
Region  2:  Bricks  rotation 
Region  3:  Two  bricks  contact 


Region  1 :  Linear  region 
Region  2:  non-linear  region 


•  Mechanisms  for  brick-and-mortar  and  laminar  composites  are  quite 
different 

•  Brick-and-mortar  Young’s  modulus  >  laminar  Young’s  modulus 

•  Elastic  properties  are  dependent  on  size,  especially  for  brick-and-mortar 
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Uniaxial  compression  and  tension  tests 


Brick  &  Mortar  Compression  Test 


Brick  &  Mortar  Tensile  Test 
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•  Pressure  along  the  compression  and  tension  axis  given 

•  Response  to  tension  is  more  linear  than  response  to  compression 

•  Off-axis  motion  of  bricks  results  in  nonlinear  response  to  the  compression 
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Surface  roughness 
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Surface  roughness  implemented  by  randomly  deleting  surface  atoms  of  the 
ceramic  layers  (polymer  atoms  are  not  shown) 
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Large  scale  molecular  dynamics  simulations  of  hybrid 
materials  with  smaller  fraction  of  polymer 


•  Large  scale  MD  simulations  of  hybrid  materials  to  understand  the  effect  of  surface 
roughness  and  functionalization  to  the  mechanical  properties. 

•  System  size  extends  to  32.5  nm  in  x  and  y  axis  and  25.7  nm  in  z  axis. 

•  Polymer/Ceramic  ratio  is  95%. 

•  Total  number  of  atoms  is  3,000,000. 
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